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Abstract. We use detailed numerical simulations and 
the v Andromedae planetary system as a template to 
evaluate the capability of the ESA Cornerstone Mission 
GAIA in detecting and measuring multiple planets around 
solar-type stars in the neighborhood of the Solar Sys- 
tem. For the outer two planets of the v Andromedae 
system, GAIA high-precision global astrometric measure- 
ments would provide estimates of the full set of orbital 
elements and masses accurate to better than 1-10%, and 
would be capable of addressing the coplanarity issue by de- 
termining the true geometry of the system with uncertain- 
ties of order of a few degrees. Finally, we discuss the gen- 
eralization to a variety of configurations of potential plan- 
etary systems in the solar neighborhood for which GAIA 
could provide accurate measurements of unique value for 
the science of extra-solar planets. 

Key words: Astrometry - Stars: individual (v An- 
dromedae ) - Planetary systems - Space vehicles - Meth- 
ods: data analysis 



1. Introduction 

As of today, five normal stars in the solar neighbor- 
hood are known to harbor either two or more planets 
(Butler et al. 1999, Mayor et al. 2000, Marcy et al. 2001, 
Udry et al. 2000) or systems composed of a planet 
and a brown dwarf (Marcy et al. 2001, Udry et al. 2000, 
Els et al. 2001). The observational evidence of the first 
extra-solar planetary systems has immediately raised cru- 
cial questions regarding their formation and evolution. Are 
the orbits coplanar? Are the configurations dynamically 
stable? The quasi-coplanarity hypothesis is one of the dog- 
mas in present planet formation theories (Lissauer 1993). 
For example, in the case of the three-planet system around 
the F8V star v Andromedas (v And hereafter), Laughlin 
& Adams (1999) and Stepinski et al. (2000) have shown 
that orbital evolution and long-term stability of the sys- 
tem can be greatly affected by small variations in the 

Send offprint requests to: A. Sozzetti (alex@phyast.pitt.edu) 



three-companion orbital model. In particular, studies on 
the dynamical stability of the outer two companions re- 
quire constraints to be placed on the maximum allowed 
values for the masses and on the range of allowed rela- 
tive inclinations of the two orbits. Radial velocity mea- 
surements cannot determine cither the inclination i of the 
orbital plane with respect to the plane of the sky or the 
position angle f2 of the line of nodes in the plane of the 
sky: without knowledge of the full set of orbital parameters 
and true mass values, general conclusions on the architec- 
ture, orbital evolution and long-term stability of the newly 
discovered planetary systems may be questionable. 

In this Letter, we gauge the potential of GAIA in de- 
tecting and measuring multiple planets in favorable cases, 

1. e. well-spaced, well-sampled orbits, and high "astromet- 
ric" S/N ratios (see Sect. 3). In order to do so, we focus on 
the v And system, which, given its characteristics, consti- 
tutes a well-suited laboratory. In this exploratory Letter, 
we do not address issues such as insufficient orbital sam- 
pling, resonant configurations, and astrometric signatures 
close to the detection limit, which clearly need to be con- 
sidered in a more realistic simulation. Work on these issues 
is in progress and will be presented in the future. 

2. Data simulation and analysis method 

The code for the reproduction of GAIA global as- 
trometric observations and analysis of the simulated 
dataset has been thoroughly described in our previous 
work (Lattanzi et al. 2000), where we defined detectabil- 
ity horizons and limits on distance for accurate orbital 
parameters and mass determination in the case of single 
Jupiter-mass planets orbiting single solar-type stars. 

We have modified our simulation code to allow for 
the presence of multiple planets around each target star 
(Casertano & Sozzetti 2000). The stellar motion is de- 
scribed in terms of the 5 basic astrometric parameters and 
of the gravitational perturbations produced by the orbit- 
ing planet(s), which are assumed to be linear (i.e., as the 
sum of independent Keplerian motions). Such simplifica- 
tion (no mutual interaction terms between planets) has 
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little impact on the significance of the analysis. In fact, 
on the time-scale of GAIA observations (5 years), vari- 
ations of the orbital parameters can be confidently con- 
sidered negligible. For planet detection we employ a stan- 
dard x 2 -Test with confidence level set to 5% applied to 
the observation residuals (Lattanzi et al. 2000), while in 
the case of multiple orbits the iterative method for the so- 
lution of the non-linear systems of equations of condition 
is based on a slightly modified version of the Levenberg- 
Marquardt algorithm (Casertano & Sozzetti 1999), which 
ensures stability of the solution. 

3. GAIA observations of the v And system 

In addressing the problem of the detection and measure- 
ment of multiple planets with GAIA, we have utilized the 
v And planetary system as a representative test case. The 
results are discussed below. 

3.1. Detection 

GAIA (for a detailed review of the mission con- 
cept see Gilmore et al. 2000) surveys the sky utiliz- 
ing a Hipparcos-like scanning law. Its elementary uni- 
dimcnsional measurement is the longitude tp of the Ref- 
erence Great Circle being scanned at a given time. The 
single-measurement astrometric error a^p was fixed to 10 
/uas throughout all our simulations: this value applies to 
stars brighter than V = 12 mag. Sensitivity to planetary 
perturbations is driven by the astrometric error cr^ of each 
measurement, which has to be compared to the magnitude 
of the gravitational perturbation induced on the parent 
star by one or more orbiting planets: the astrometric sig- 
nature a (in arcsec) is the apparent amplitude of the or- 
bital motion of the star around the system barycenter, 
defined as 



MpOp 
M* d ' 



(1) 



where M p , M+ are the masses of the planet and star re- 
spectively (in solar masses), a p the semi- major axis of the 
planetary orbit (in AU), d the distance of the system from 
the observer (in pc). 

We simulate 200 planetary systems on the sphere us- 
ing the parameters of the v And system in the original 
paper by Butler et al. (1999), as summarized in Table 1. 
Two orbital elements are not determined by radial velocity 
measurements, namely i and Q. We simulate systems with 
random uniform distribution in O and express detection 
probabilities as a function of the inclination of the orbital 
plane (5° < i < 90° ). 1 We assume the orbits of the three 
planets are perfectly coplanar, i.e. i and SI are exactly the 
same for all the system components. Planet masses are 
scaled with the co-secant of the inclination angle. During 

1 Astrometric detection probability does not depend on the 
apparent orientation of the orbit (Lattanzi et al. 2000) 



Table 1. v Andromeda: stellar parameters and orbital ele- 
ments utilized in the simulations. Planet masses and astromet- 
ric signatures are computed as lower limits corresponding to 
edge-on configurations (sini = 1) 



Parameters 



v And 



A (deg) 




37.9 


13 (deg) 




28.9 


(mas/yr) 




-242.74 


na fmas/vr) 




-266.98 


n (mas) 




74.25 


Visual Magnitude 




4.09 


Spectral Type 




F8V 


d (pc) 




13.47 


M* (M ) 




1.3 


Orbital elements v And B 


v And C v And D 


q (/uas) 2.40 


100.1 


658.4 


a (AU) 0.059 


0.83 


2.50 


P (d) 4.6 


241.2 


1266.6 


e 0.034 


0.18 


0.41 


t (JD) 24500088.6 


2450154. 


9 2451308.7 


uj (deg) 83.0 


243.0 


247.7 


MpSini(Mj) 0.71 


2.11 


4.61 



the 5-year mission lifetime, the system is observed for a 
total of 70 times 2 . 

We find that, given the favorable values of the "as- 
trometric" signal-to-noise ratio S/N — a/a^, the pres- 
ence of both the outer and middle planet can be easily 
detected (at the 95% confidence level) by GAIA (in the 
residuals to a single-star model and to a model contain- 
ing the outer companion, respectively), regardless of the 
inclination of the orbital plane. In particular, the system 
configuration is such (orbits with well sampled periods, 
high S/N, no resonances) that detectability of the middle 
planet is not affected by the presence of the outermost 
companion. On the other hand, the inner planet, v And 
B, cannot be detected, unless i < 10°: as a matter of fact, 
its orbital period, shorter than the average observing pe- 
riod, is poorly sampled by GAIA, and most importantly 
the magnitude of the astrometric signature is such that, 
for i > 10°, S/N < 1 and detection through the x 2 -Tcst 
fails. 

3.2. Orbit reconstruction 

We focus here on GAIA's ability to make accurate mea- 
surements of the orbital elements and masses for the outer 
two planets of the system, which produce the greatest as- 
trometric perturbations. Specifically, we have simulated 
systems containing the middle or the outer planet only 
(utilizing a single-planet solution for the orbital fits), and 



2 The GAIA scanning law is such that the number of basic 
observations is function of the ecliptic latitude [3 
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the complete configuration with the three planets, fit- 
ting to the observation residuals a) a two-planet model, 
disregarding the presence of the gravitational perturba- 
tion induced by the inner planet, or b) the complete 
three-planet model (for details on the fitting procedure, 
see Casertano & Sozzetti 1999 or Sozzetti et al. 2001). 
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Fig. 1. RMS errors for M p and i, in the case of the outer 
two planets of the v And system, expressed as fraction (%) of 
the true values of the parameters: for v And D, mass can be 
measured to 0.1-1% accuracy and inclination to better than 
1% (for i > 30°); for v And C, mass is accurate to 1-3% and 
inclination to better than 10% (for i > 30°) 



In Figure 1 we have summarized the results for two of 
the most important parameters of the system, inclination i 
(directly measured) and planet mass (computed once a re- 
liable estimate for the mass of the parent star is provided), 
for v And C and v And D. Owing to the difference in as- 
trometric signature (see Table 1), the outermost planet 
is measured more accurately than the inner planet (ac- 
curacy is almost an order of magnitude higher). 3 As the 
inclination i decreases, towards a face-on configuration, 
the true mass inferred for each planet from the radial ve- 
locity data increases, therefore the astrometric signature 
is larger and provides a more accurate measurement of the 
planet's mass. At the same time, the observations are less 
sensitive to the inclination itself, and its fractional mea- 
surement error increases. It should also be noted that, if 
the system is very close to face-on, the accuracy of the in- 
clination measurement will be substantially increased by 
combining radial velocity and astrometric data. 



Finally, we emphasize the fact that, similarly to what 
happened for detection, the orbital arrangement of the v 
And system is such that it does not affect our ability to 
reconstruct each planet's orbit. As shown in Figure 1, the 
two outer planets can have orbital elements and masses 
measured almost as well as if each of them was the only 
orbiting companion. For a two-planet orbital fit to a com- 
plete three-planet simulation, the quality of orbit recon- 
struction degrades only for i < 10°, when the signal from 
v And B becomes greater than the single-measurement 
error (we recall = 10 /ias) , and thus detectable in the 
observation residuals. 



3.3. Coplanarity analysis 

As discussed in the Introduction, the unexpected orbital 
configuration of the v And system has motivated detailed 
theoretical studies on its dynamical evolution and long- 
term stability. For example, neglecting to first order the 
effects of the innermost planet on the overall stability of 
the system, Stepinski et al. (2000) come to the conclusion 
that dynamical stability requires the orbital inclination 
of the outer two companions to be greater than i ~ 13°, 
otherwise the two objects would be too massive and grav- 
itational interactions would disrupt the system. Further- 
more, the system cannot be stable in the long term if 
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3 Note that the multiple-planet fit provides very accurate 
estimates for all of the orbital elements of the outer two planets. 



20 40 60 80 100 

Orbital Inclination (deg) 

Fig. 2. Coplanarity analysis for the v And system: for i > 20°, 
the orbits of the two outer planets can be confidently identified 
as quasi-coplanar, with (absolute) uncertainties of order of a 
few degrees (errors are computed using the formal expressions 
from the covariance matrix of the multiple-planet fit) 
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relative inclinations are greater than 55°, 35°, and 10°, 
for i ~ 64°, i ~ 30°, and i ~ 15°, respectively: then, the 
more massive the planets, the closer to coplanarity their 
orbits have to be, for the system not to be destabilized on 
a short time-scale. 

Figure 2 shows the estimated accuracy with which 
GAIA would measure the coplanarity of the orbits of v 
And C and v And D (as shown in Sect. 3.1, unless its orbit 
is almost face-on, v And B is not even detectable). Fitting 
a three-planet model to the simulated observations, we es- 
timate differences between position angles of the lines of 
nodes (A^cd) < 10° and differences between inclinations 
of the orbital planes (A«cd) < 3°, for i > 20°. Thus, quasi- 
coplanarity can be reliably established. As i decreases, the 
uncertainty on Af^cD grows, because of the increasing dif- 
ficulty in the correct identification of the position angle f2, 
which becomes undefined at i = 0°. Simulations of non- 
coplanar systems with AficD ^ 0° and/or Aiqd ^ 0°, 
yielded similar results. Our findings indicate that GAIA 
high-precision measurements would be instrumental in 
verifying the system stability, except for quasi-face on 
orbits: however, the inclination of v And D is probably 
greater than 10°-15°, otherwise the astrometry of Hippar- 
cos would have revealed its presence (Mazeh et al. 1999). 

4. Discussion and conclusions 

The ability to accurately measure orbital parameters and 
masses for a potentially large number of planetary sys- 
tems with massive planets is the key element of the GAIA 
contribution to the science of extra-solar planets. 

Utilizing the v And three-planet system as a tem- 
plate, we have shown how, in case of planetary sys- 
tems with planets producing favorable astrometric signa- 
tures (S/N > 10) and with well-sampled orbital periods 
(shorter than the 5-year mission lifetime), the parameters 
and masses of each companion can be measured almost as 
well as if it was the only planet. 

We have made coplanarity tests, and found that GAIA 
would be able to determine whether orbits are coplanar 
or not with uncertainties of a few degrees, thus providing 
theory with the observational evidence needed to address 
the long-term stability issue for the system . 

Lattanzi et al. (2000) derive a 60-pc distance horizon 
for 30% accuracy measurements of orbital elements and 
masses in the case of single Jupiter-mass planets orbiting 
solar-type stars with periods ranging between 0.5 and 11.8 
years (the true period of Jupiter). 

We can generalize those findings, in an early attempt 
to identify the variety of possible configurations of plane- 
tary systems detectable and measurable with GAIA. The 
results reported in this letter allow us to say that the 
60-pc limit on distance holds for detection and measure- 
ment accurate to 30%, or better, of planetary systems com- 
posed of planets with well-sampled periods (P < 5 yr), 
and with the smaller component producing an astromet- 



ric S/N > 2. Accurate coplanarity tests will be possible 
for systems producing S/N > 10. 

The frequency of multiple- planet systems, and their 
preferred orbital spacing and geometry are not currently 
known. Based on star counts in the vicinity of the Sun ex- 
trapolated from modern models of stellar population syn- 
thesis, constrained to bright magnitudes (V < 13 mag) 
and solar spectral types (earlier than K5), we should 
expect <~ 13 000 stars to 60 pc (Lattanzi et al. 1999). 
GAIA, in its high-precision astrometric survey of the so- 
lar neighborhood, will observe each of them, searching for 
planetary systems composed of massive planets in a wide 
range of possible orbits, making accurate measurements of 
their orbital elements and masses, and establishing quasi- 
coplanarity (or non-coplanarity) for detected systems with 
favorable configurations. The size of the stellar sample is 
such to ensure that the results obtained by GAIA would 
constitute a fundamental complement to those which will 
come from other planet-search programmes, in order to 
build the necessary statistics for deeper theoretical under- 
standing of planetary formation processes. 
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